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Abstract 



G I The radiative B -^ t^t^^ decay is investigated in the framework of the two Higgs 

doublet model . The dependence of the differential branching ratio on the photon energy 

and the branching ratio on the two Higgs doublet model parameters, mff± and tan/3, are 

CSI . studied. It is shown that there is an enhancement in the predictions of the two Higgs 

doublet model compared to the Standard model case. We also observe that contributions 

1^ • of neutral Higgs bosons to the decay are sizable when tan (3 is large. 
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1 Introduction 

Rare B-meson decays are one of the important research areas to test the theoretical models 
and make estimations about their free parameters. In the Standard model (SM) they are 
induced by flavour changing neutral currents (FCNC) at the loop level. This ensures a precise 
determination of the fundamental parameters of the SM, such as Cabbibo-Kabayashi-Maskawa 
(CKM) matrix elements, leptonic decay constants, etc. In addition, the studies on rare B- 
meson decays give powerful! clues about the existence of model beyond the SM, such as two 
Higgs doublet model (2HDM), minimal supersymmetric extension of the SM (MSSM) |^, etc. 
Among rare B-decays, B -^ i^i~'y decays are of special interest due to their cleanliness and 
sensitivity to the new physics. They have been investigated in the framework of the SM 
in 0, Q for £ = e, /i and in [Q for £ = r. The theoretical results given in ^ and [Q are 
BR(5, -^ e+e-7) = 2.35 X 10-9, BR(5, -^ AtV"7) = 1-9 x lO^^ and BR(5, -^ r+r-7) = 
9.54 X 10~9, respectively. These decays get negligible contributions from the diagrams, where 
photon is radiated from any charged internal line due to the fact that they will have a factor 
ml/M^^ in the Wilson coefficients. When photon is radiated from the final charged leptons, the 
contribution is proportional to the lepton mass rrii. Therefore, for £ = e, /i case, it is negligible; 
however for £ = r it gives a considerable contribution to the amplitude. In the 2HDM, there is 
a part coming from exchanging neutral Higgs bosons and in contrast to i? ^ i^i^'j {£ = e, /i) 
decays , we could expect that they significantly contribute for Bg -^ t^t^'j decays. Therefore, 
in this work we study the B^ — > r^r^^ process in the framework of the 2HDM (Model I and 
II). 

2IIDM is one of the simplest extensions of the SM, obtained by the addition of a second 
Higgs doublet. In this model, there are one physical charged Higgs scalar, two neutral Higgs 
scalars and one neutral Higgs pseudoscalar. The Yukawa lagrangian causes that the model 
possesses tree-level FC couplings of the neutral Higgs particles. To avoid such terms, it is 
proposed an ad hoc discrete symmetry on the 2HDM potential and the Yukawa interaction. 
As a result, it appears two different choices for how to couple the quarks to the two Higgs 
doublets: In the first choice (Model I), the quarks do not couple to the first Higgs doublet, but 
couple to the second one. In the second choice, (Model II), the first Higgs doublet couples only 
to down -type quarks and the second one to only up-type quarks. 

The paper is organized as follows: In sec. 2, we present the theoretical framework for the 
Bg — > r"'"r~7 decay and describe some details of its decay rate calculation . We give a numerical 
analysis and discussion of our results in sec. 3. Appendices contain a list of the operators and 
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the Wilson coefficients, as well as some relevant formula about the long distance contributions. 

2 Bs ^ T^r~j decay in the framework of the 2HDM 

The exclusive decay Bs —>■ r^r^'-f can be obtained from the inclusive one b —>■ st^t^'j. In 
order to calculate the relevant physical quantities for the decay b -^ st^t"'-/, we start with the 
QCD corrected amplitude for the process b -^ st^t^ . At this stage, the effective Hamiltonian 
is obtained by matching the full theory with the effective low energy one at the high scale 
yU. The Wilson coefficients are evaluated from /i down to the lower scale fi ~ 0{mh) using the 
renormalization group equation (RGE). The effective Hamiltonian in the 2HDM for the process 
b -^ ST^T^ is 

In this equation Oj are current-current (z = 1,2), penguin (i = 1,..,6), magnetic penguin 
(i = 7, 8) and semileptonic (i = 9, 10) operators . The additional operators Qi, {i = 1, .., 10) are 
due to the neutral Higgs boson exchange diagrams, which give considerable contributions in the 
case that the lepton pair is t^t^ 0. Cj(/i) and Cq^{^) are Wilson coefficients renormalized at 
the scale /x. All these operators and the Wilson coefficients, together with their initial values 
calculated at /x = mw in the SM and also the additional coefficients coming from the new Higgs 
scalars are presented in Appendices A and B. The QCD corrected amplitude for the inclusive 
b -^ sT+r- decay in the 2HDM (Model I or II) is 

M = ^V,,V:,i^Clf^{-s^,PLb)f^,r + C,o{si,PLb)fi,i,r 

- 2C'j^{sia^i,puPRb)f-i^T + CQ^{s-i^PRb)fT + CQ^{s-i^PRb)f-i5T\ . (2) 

where Pl,r = (1 =F75)/2 , p is the momentum transfer and Vi/s are the corresponding elements 
of the CKM matrix. 

In order to obtain the matrix element for b -^ st^t^'j decay, a photon line should be at- 
tached to any charged internal or external line. The contributions coming from the attachement 
of photon to any internal line are suppressed and we neglect them in the following analysis. 
We now start with the case in which a photon is attached to the initial quark lines. The 
corresponding matrix element for the Bg — > r^r^'-f decay is 

M, = {i\M\B) = ^^Vt,Vzlc',"f^,T{^\s^,{l - 75)&|5) + Cior7^75r(7|s-7^(l - l,)b\B) 



2C7^(7|3za^.P.(l + l,)b\B)f-f^r + CQ,fr{-f\s{l + j,)b\B) + CQ,f-f,r{-f\s{l + j,)b\B) 



(3) 



These matrix elements can be written in terms of the two independent, gauge invariant, 
parity conserving and parity violating form factors p, 0: 



(7|s7/.(lT75)&|5) = —T{^f^apae*aPpqa9{p'^)±iel{pq)-{e*p)q^ f{p 



m 



B 



(4) 



and 



{-i\sia^,^p^{lT'lb)h\B) 



TTl, 



e^^apaelPisqa giip"^) T i C{pq) - {^*p)q^^ h{p 



(5) 



Here e^ and q^ are the four vector polarization and four momentum of the photon, respectively 
. To calculate the matrix elements (7|s(l±75)6|-B), we multiply both sides of eq. (||) by p^ and 
use the equations of motion. However, neglecting the mass of the strange quark they vanish. 



(7|s(l±75)&|5) = 



(6) 



Substituting Eqs. (^) and (|^) in (^, for the matrix element A^i (structure dependent part) we 
get 

aGp ' 



Ml 



2^2 



vr 



-VibV^ei ef,at3aeaPf3qa [A r-f^r + C r7^75r] + i e*{pq) - {e*p)q^ [Br^i^T + DT-f^-f^T 



(7) 



where 
A 



ml 



C^," gip') - 2C,^ giip 



pz. 



B 



m 



B 



yeff 



nib 



pZ 



C = ^9{p') , /^ = %/(/) . (8) 

m% m^ 

Note that the neutral Higgs exchange interactions do not give any contribution when pho- 
ton is attached to the either one of the initial quark lines. However, when a photon is radi- 
ated from the final r-leptons the situation is different and the corresponding matrix element 
(Bremsstrahlung part) is 



M, 



+ 



aGf 

2v^7r 

rUb 



VtbV*etfB{ 2m,Cio + ^ia 



2 ?t1t 



1 1 

+ 



2piq 2p2g, 
3 



nib J [ \2piq 2p2qJ 

T fT + T \ r 

\2piq 2p2qJ 



I5T 



(9) 



where we have used 



(0|s7^75fe|5) 
{0\sa^,{l±^,)b\B) 

(0|s75&|5) 
{0\sb\B) 



ifsP. 



Bfi , 



. 



rrii 



nib 



(10) 



and the conservation of the vector current. Here Pb is the momentum of the B-meson. 
Finally, we get the total matrix element for the B -^ t^t^'-/ decay as 

M = M1+M2 ■ 



'11^ 



To calculate the decay rate, we need the square of this matrix element. By summing over the 
spins of the r-leptons and the polarization of the photon, we obtain 



lA^I 



\Mi\^ + \M2\^ + 2Re{MiMl) 



'121 



where 



\Mi\ 



2Re{MiM*2] 



aGi 



2v/2 



TT 



-Vti>v:, 



+ 4 \C\^ + \D\^ 



+ 4 



\A\ 



\B 



47ra <^ 8 Re {B*C + A* D) / {piq - p2q) {piQ + P2q) 
{piqf + {p2qf) - 4m^ (piq) {p2q) 
2^ ^{piqf + {P2qf) + 4m^ (piq) (pa?) 



p^ — 2m^ 



p 



2m: 



aGt 



2^2 



vr 



-VttVZ 



Anal leCio/sm^ 



Re (A) 



{Piq + P2qf 
(Piq) {P2q) 



+ Re{D) 



{piq + P2qf {piq - P2q) 



ij>iq) ij>2q) 



m. 



rub 



m. 



■Cqi 



Re{B) 



{piq + P2qY 
(Piq) {P2q) 



- Re(C) 



{piq + P2qf' {piq - P2q) 



+ -^ReiB) 
rrib 



^{ ml -3p2q){piq) {2ml - P^){P2q) 



P2q 



Piq 



(piq) ij>2q) 



(13) 



(14) 



\M, 



-16 



aGi 



2V2 7r 
2mrCio 



VnV* 



tbVts 



Ana 



Uli 



■Cr 



rrib 



+ 



2mi 



{piqY 



m 



B. 



'^ICqi 



2mrCio + ^^Cq, 
^b 



rrib 



p^ + 2p2q 



'm\CQ^' 



rrib 



p^ + 2^20' — Ami 



+ 



+ 



4 
piq 

2ml 



m 



2mrC,o + ^^CQ, + 



'm|C, 



TUb 



nib 



+ 



{P2q) 
4 
P2q 



2mrCw + ^CQA (p^ + 2p,q 



3m^ — p — 2p2q 

(^^^^^]\p' + 2p,q-4ml 
\ rUb J ^ 



^Br^ \ , (^%Cqi 



2m,C,o + ^^CQ, + 

Vflb 



nib 



3nil — p^ ~ 2]9ig 



{Piq) {P2q) 



2 \ ^ / 2 /^ 

^^B r^ \ 2 /o™2 „2\ P^B^Qi 



2ni^Cio + ^^Cq, p' {2nii - p') - 

\ ^6 / 



TUb 



1 (p'^ + 2p2q - Aml^ ] } . 

(15) 



Here pi, p2 are momenta of the final r-leptons. 

In the rest frame of the i?-meson, the photon energy E-y and the lepton energy Ei are 
restricted in the region given by 

ml, — Ami 

0<E^< -^ ^ , 

2mB 



niB — E^ 



E^ 



2\ 



1- 



Aml 



E^ 



m| - 2mBE^ - i - ^ 2 \ 



1 



Ami 



m\ — 2mBE^ 



(16) 



In \M.2\ it appears an infrared divergence, which originates in the Bremstrahlung processes 
when photon is soft and in this limit, the B -^ r^r^'-f decay cannot be distinguished from 
B -^ T^T^ . Therefore, in order to cancel the infrared divergences in the decay rate both 
processes must be considered together. In ref. [^ it has been shown that infrared singular 
terms in \M.2\ exactly cancel the 0{a) virtual correction in S -^ t^t~' amplitude. However, 
in this work we consider the photon m B ^ t^t^'-/ as a hard photon, following the approach 
described in ref. [^ and impose a cut on the photon energy. The lower limit of this cut is 
choosen so that the radiated photon can be detected in the experiments, namely E^ > 50 MeV 
(^ a niB with a > 0.01). After integrating over the phase space and taking into account the 
cut for the photon energy we get for the decay rate 

2 



aGf 



2^2 



TT 



-VtbV,: 



a 



(2vr) 



ni^TT 



m 



B_ 

12 



l-4r 



x-^dx \ 1 



Ar 
1 — X 



4/ijr{ ( Cio + 



2mbmr ^ 



l-4r 



{\Af + \Bf) (1 - a; + 2r) + (\Cf + \Df) {1 - x - 4r) 
l + i/l 



4r 



x^ dx Re {A) In- 



1 — a; 



4r 
1 — X 



+ 



m, 



l-4r 



^—Cn, I xdxRe{B) 

2mhmr Js 



:i-x) wi 



4r 
1 — X 



'1- 



+ (x - 2r) 



4r 
1 — X 



1-A 1 



4r ■ 
1 — X 



mi, ^ \ Znihrn^ 



Q2 



l-4r 



dx 



[1 — x) / 4r 



x 



1 — a; 



,, 2r 1 \ +Y^~r^^ 

+ (IH X In V 

X X / / 4r 



1-A 1 



1 — X 



1 / mBCq ,^^ /■1-4'- 
r V 2mb 



(ix 



, (1 - x) / 4r 



1 — X 



4r 



.2 1 



-iH \ hx + -(4x - 6) In V 

XX X / / 4r 



1 — X 



17) 



where r = — ^, 5 = 2a and x = — - 



m 



B 



rriB 



is the dimensionless photon energy satisfying 



S <x<l 



Arriz 



mi 



In our numerical calculations, we use the dipole forms of the form factors given by 



9(P^) 



IGeV 



„2 

P \2 



fiP') 



0.8 GeV 



5.62 



(1 - ^f 



P ^2 
6.52' 



9i{p^ 



3.74 GeV^ 



W) 



40.5' 



0.68 GeV^ 
^ 30^ 



;i8) 



which were calculated in the framework of the light-cone QCD sum rules [^ ^. 

3 Results and Discussion 

In the 2HDM there are number of free parameters, namely masses of the charged and neutral 
Higgs bosons {m,fj±, m^^o, m^o), the ratio of vacuum expectation values of Higgs bosons, tan/5 = 
f2/fi, and the angle a due to the mixing of neutral Higgs bosons, A^ and h^. The values of 
these parameters are restricted by using the existing experimental data. The non-observation 
of charged if ^ pair in Z decays [§ gives the model independent lower bound of the mass of 
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the charged Higgs H^, f^tR^ ^ 44 GeV. However there is no experimental upper bound for 
mH± except mH± < 1 TeV coming from the unitarity condition |]10[. Further, top decays 



give mff± > 147 GeV for large tan/? |jTT[. The other parameter of 2HDM, tan/?, is restricted 



as tan (3 > 0.7 from Z ^> hh decay |12|. The ratio tan/^/m/zi can also be restricted and 



it has been estimated as tan/5/r7i/f± < 0.38 GeV^^ |T^ and icml3/mH± < 0.46 GeV^^ |]T4 
from the experimental results of the branching ratios of the decays B ^ tv and B ^ X tv. 
The upper bound has also been given for the same ratio as ta,'n.(3/mH± = 0.06 GeV~^in the 
case that sufficient data could be taken and the theoretical uncertainties could be reduced for 



the exclusive decay B ^ Dtu [^ . Recently, the relation between mH± and tan /? has been 



estimated in [|T6[ , taking into account the CLEO measurement of the decay B -^ Xg^y [0 



Br{B -^ X,7) = (3.15 ± 0.35 ± 0.32) 10"* , (19) 

In our calculations, we take the masses r/i/iO and m^o equal and not too heavy since the b- 
quark dipole moment is strongly sensitive to the difference between these masses in the 2HDM 
|T8| . Further, we choose the value of the angle a as being zero since the mixing between h^ and 
A^ is weak. For completeness, we have also checked the dependence of the branching ratio on 
a for the fixed values the other 2HDM parameters and seen that this dependence is negligible. 

In the present work, we study the 2HDM parameters dependence of the BR and dimension- 
less photon energy dependence of the differential branching ratio (dBR/dx) in Model I and II. 
Doing this, we have used the input parameters given in Table I. 

In Fig.l, we present dBR{B -^ T^T^'~f)/dx as a function of x = 2E^/mB in the SM and 
in Model II for niHi = 400 GeV and tan/? = 2. We do not display the predictions of Model 
I there, since they are very close to those of Model II. In this figure, curves with sharp peaks 
represent the long distance contributions. From Fig. 1, we see that there is an enhancement in 
the 2HDM compared to the SM case. 

Fig 2. shows the dependence of the BR on the Higgs boson mass mH± for different values of 
the parameter tan/? for Model I and II, as well as for the SM. We again observe an enhancement 
for the BR in 2HDM compared to the SM case. For example, for niH^ = 400 GeV and tan/3 = 2 
, BR{B -^ r+r-7) = 4.18 X 10^^ in Model I, and BR{B -^ r+r-7) = 4.20 x 10'^ in Model H. 
These values are greater than the SM predictions, which is BR{B — > r^r~7) = 4.13 x 10~®. In 
addition, the mH± dependence of the BR becomes weaker with increasing values of tan/? for 
both models. 

We present the BR as a function of tan (3 for different values of mH± in Model I and II in 
Fig. 3 and Fig. 4, respectively. It is seen that additional contributions coming from neutral 



Higgs exchange diagrams (i.e., contributions with Cq. ^ ) causes the BR to increase with 
the increasing values of tan/3 in contrast to the case that neutral Higgs do not contribute 
{Cq. = 0). The reason for these two different behaviors can easily be understood by comparing 
eqs.(^) and (|22|) , which represent the neutral Higgs bosons and the remaining contributions, 
respectively, namely the first one is proportional to tan^/? , while the second is l/tan^/5 so 
that for the larger values of tan /3 , neutral Higgs contributions dominate in the BR. 

As a conclusion, we observe an enhancement in the differential branching ratio and the 
branching ratio of the exclusive process B -^ r^r^'y in the framework of the 2HDM as compared 
to the SM. Further, this enhancement becomes more detectable for large tan/? values lying in 
experimentally restricted regions. Therefore, the measurement of this exclusive decay gives 
important clues about the new physics beyond the SM, corresponding model parameters and 
also the effects of neutral Higgs contributions. 
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A The operator basis 



The operator basis in the 2HDM (Model I and II) for the process under consideration is [0, ^ 



O2 = {sLa1^lCLa){cLpl^hL|3), 

q=u,d,s,c,b 

O4 = {sLa^tiblls) Yl iQL/3YqLa), 
q=u,d,s,c,b 

q=u,d,s,c,b 

Oe = (sLalf^bLp) Yl {QRHI^QRo)^ 

q=u,d,s,c,b 



e 



■fiu 



g _ 

'^S = T7r^{sLal,j,bLa){lYl) , 
g _ 

0^0 = -^^^{sLalphLa){ll^lbl) , 



afiu 



e 



^1 = T7r^(^lb%) (fr) 



IQtv 



Q2 = 77^(^1 &S)(^75r) 



167r^ 
Qs = ^M^r) E (IlQ^r) 

q=u,a,s,c,b 
q=u,d,s,c,b 
g=M,d,s,c,fe 

Qe = ^.(-^I^r) E (?"^^2) 

q=u,d,s,c,b 
q=u,d,s,c,b 



Q9 = Yf^(^l^''^4) E {fL<y,^<fR) 

q=u,d,s,c,b 

Qio = ^2(''l^"'bR) E i^'R^,^Q2) (20) 

q=u,d,s,c,b 

where a and (3 are St/ (3) colour indices and JF^*^ and ^'^'^ are the field strength tensors of the 
electromagnetic and strong interactions, respectively. 

B The Initial values of the Wilson coefficients. 

The initial values of the Wilson coefficients for the relevant process in the SM are [pi 

Cl,3,.. . 6,11,12 ("^w/) = , 



C^'imw) = 1, 

— -^/ CAT"/ \ Oi/y ^Jj (Ji/y ^dU ~I~ I lAj 

Cfo"(n>»-) = -^TT (-Bl^) - C(i)) , 
sm 6w 

C^f'imw) = ^ = 1,..,10. (21) 

The initial values for the additional part due to charged Higgs bosons are 

cZ.A^w) = 0, 

C^imw) = XFM +YFM , 

Ciimw) = XGM +YGM. 

C^{mw) = XH,{y) , 

C^,{mw) = XL,{y) , (22) 

and due to the neutral Higgs bosons are P 

+ (sin^ a + hcos^ a)fi{x,y) + [m^o/m^ + (sin^a + hcos^ a){l — z)]f2{x,y) \ 

10 



Cgimw) 



rribe 






rrirg^ 

,i = 5,..,10 , 



(23) 



where 



m: 



mr 



X 



X 



m 



y 



w 



m 



H 



h 



m 



ho 



and 



f2{x,y) 



X 



xlny 



y 

Inz 



{z-x){x-l) {z-l){x-l) 

1 / 1 \ y -1 



m 

, My) 



, fiix,y) 



xlnx ylny 



HO x-1 y~l ' 

l-y + y\ny 



iy-i] 



(24) 



(1) in Model I (//) 



Fi{y) 
F2{y) 
Gi{y) 
G2{y) 
Hi{y) 
Li{y) 



tan^/3 \ta.\i^ 13 j ' tan^/5 

The exphcit forms of the functions Fi(2)(i/), G'i(2)(y), Hi{y) and Li{y) are given as 

y{7-5y-8y^) , y\3y - 2) 



(25) 



72(y-i)3 +i2(y-i)4 ""y^ 

y{5y-3) y{-3y + 2) 
12(y-l)2+ 6{y-ir ''^' 



|/(-|/^ + 5|/ + 2) -y^ 

24(t/ - 1)3 A{y - 1)4 

y{y-3) y 

4{y - 1)2 2{y - 1 
1 — Asin'^Ow xy 



;^lni/, 



sin^d 
1 



H^ 



Iny , 



1 



xy 



sin?6w 8 



y-l (l/-l)2 

1 



Iny 



47|/2 - 79?/ + 38 3|/3 - 6?/ + 4 



108(1/ -1)3 



18(1/ -1)4 



Iny 



+ 



y-l (y-1^ 



Iny 



(26) 



Finally, the initial values of the coefficients in the 2HDM are 



a 



2HDM 



(mw) = C^''{mw) + C^{mw) 



(27) 



Using these initial values, we can calculate the coefficients Cf^^'^{fi) and CQ^^^\fi) at any 
lower scale in the effective theory with five quarks, namely u, c, d, s, b similar to the SM case. 
Wilson coefficients C2^^^(/i), C2^^*^(/i),C2o^^^(^), C2f^^(/^) and C2f«*^(/^) play the es- 
sential role in this process and the others enter into expressions due to operator mixing. For 
completeness we would like to give the explicit expressions of the cofficients essential in this 
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process. The effective coefficient Cj^{jj) is defined as ||2^ 



,Tnr ^Of/DAT/ X ,. ^^c ^2HDM I 






(28) 



where the leading order QCD corrected Wilson coefficients Cj 



LO,2HDM 



(/i) are given by |]T^, 



23 



^7 



LO,2HDM 



(/^) 



^16/23^2/fDAf ^^^) ^ (8/3)(r/i^/23 - r]''/'')C',''''^'{mw) 



+ Cl 



2HDM , 



mw) J2 ^iV"' 



(29) 



i=l 



and 1] = as{'m^r)/as{fJ'), hi and Oj are the numbers which appear during the evaluation [03 
The perturbative part of the Wilson coefficient Cg (/i) can be defined as |^, ^ 



crif^) 



r^2HDM 



{fi)fj{s) 



+ hiz, s) i3C^i^l) + C2ifi) + 3Csifx) + C,{fi) + ^M + CM) 

- hi{l, s) {ACM + 4C4(/i) + 3CM) + CM) 

- ^/i(0, s) (C3(^) + 3C4(^)) + ^ (3C3(/i) + C4(/i) + 3C5(/i) + C6(/x)) . 



(30) 



Here the contributions of the coefficients Ci(/i), ...., C^lfi) are due to the operator mixing. In 



eq. (30) fj{s) represents the one gluon correction to the matrix element Og with nis = pO| and 
the function h{z, s) arises from the one loop contributions of the four quark operators Oi, ...,0q 
. Their explicit expressions are 



r]{s) = l-\ iu[s) 



IT 



(31) 



where 



and 



uj(s) 



2^2 _ lLi2(s) - -lnSln(l - s) - f ^ ^^. 1n(l - s) 



9 3 ^ ' 3 

2s(l + s)(l-2s) 



3(1 



:i + 2s) 



Ins + 



3(1 + 2s) 
5 + 9s-6s2 
6(l-s)(l + 2s) ' 



(32) 



h{z, s) 



h{0,s) 



8 , mfe 8 , 8 4 
- In In z H 1 — X 

9 u 9 27 9 



9 



(2 + a;)|l-x|^/2 



In 



rr^+l 



/1-x-l 



-iir] , for X = ^ < 1 



2 arctan ■ 



-^=^, for a; = ^ > 1, 

/x— 1' s ' 



8 8, mfe 4 4. 

m m s H — tir , 

27 9 /i 9 9 



(33) 



(34) 
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where z = ^^ and s = ^. In addition to the perturbative part, there exist also the long 
distance (LD) one due to the conversion of the real cc into the lepton pair t^t~ , described by 
the reaction B -^ '-/ipi -^ 'yT~^T~ , where i = 1, .., 6. Adding this contribution to the perturbative 
one coming from the cc loop, the NLO QCD corrected Cg (/i) can be written as: 

Clf^{fl) = Cr'\fi)+Yreson{S), (35) 

where Yreson{S) in NDR scheme is defined as 

(3Ci(/i) + C2{^^) + 3C3(/i) + C,{^l) + 3C5(/i) + CM) . (36) 

The phenomenological parameter k in eq. (|36D is taken as 2.3 |24|. 
Finally, the Wilson coefficients Cq^^jj) and Cq^{h) are given by 

CQX^^) = r'^'^' CQXmw) , ^ = 1,2. (37) 
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Parameter 


Value 


Trie 


1.4 (GeV) 


rUb 


4.8 (GeV) 


^em 


137 


\VtbV* 


0.045 


rriB, 


5.28 (GeV) 


r{Bs) 


1.64 X 10-12 (s) 


rrit 


176 (GeV) 


mw 


80 (GeV) 


Vflz 


91.19 (GeV) 


m-r 


1.78 (GeV) 


TUhO 


80 (GeV) 


mHO 


150 (GeV) 


niAo 


80 (GeV) 


^^ 


mi, 


^QCD 


0.225 (GeV) 


otsijnz) 


0.117 


sinOw 


0.2325 



Table 1: The values of the input parameters used in the numerical calculations. 
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I 
+ 

t 



3M with LD 

Model II 

Model II with LD 




Figure 1: Differential branching ratio as a function of x = 2E^/mB in the SM and Model II 
for mH± = 400 GeV and tan/? = 2. In this figure, curves with sharp peaks represent the long 
distance contributions. 
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Figure 2: Branching ratio as a function of m//± in the SM, Model I and II for different values 
of tan/9. 
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Figure 3: Branching ratio as a function of tan/5 in the SM and Model I for different values of 
mH±- Curves with Cg. ^ (Cg- = ) represent the contributions including (not including) 
the neutral Higgs boson interactions. 
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Figure 4: The same as Fig |[ but for Model II. 
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